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We report on a simple model of spatial extend anti-tumor system with a fluctuation 
in growth rate, which can undergo a nonequilibrium phase transition. Three states 
as excited, sub-excited and non-excited states of a tumor are defined to describe its 
growth. The multiplicative noise is found to be double-face: The positive effect on 
a non-excited tumor and the negative effect on an excited tumor. 
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In the past decades, many studies have focused on biodynamics [1-4], specially noise 
biodynamics [5-9]. More than ever, cancer research is now an interdisciplinary effort which 
requires a basic knowledge of commonly used terms, facts, issues, and concepts. Phase 
transition of tumor growth induced by noises is one of the most novel foundations in recent 
years [10, 11]. However, in all these studies the systems are zero-dimension and insufficient 
to describe the real progress in the field of tumor growth, furthermore at present the space 
has become a fundamental variable to study [1, 12, 13]. 

Chemotherapy and Immunotherapy remain far from good understanding, although they 
as a potential practical partnership have attracted numerous attentions of scientists for at 
least one decade [14, 15]. Due to the different responses of tumor cells to chemotherapy and 
immunotherapy, more recently Lake and Robinson suggested that there is an interesting and 
significative case for combining chemotherapy and immunotherapy in tumor treatment [14]. 

In this paper, chemotherapy and immunotherapy are joined by a spatial extend anti- 
tumor model with three elements, which are (1) a spatiotemporal fluctuation of growth rate 
induced by chemotherapy, (2) an immune form, and (3) a spatial extend form. Based on the 
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analysis on its unique stochastic differential equation and relevant Fokker-Planck equation, 
we will show that the spatiotemporal fluctuation can lead to a transition of tumor-growth 
state through both theoretical analysis and numerical calculation. Although noise-induced 
phase transition is a well known phenomenon, double-faces effect of a noise on a tumor 
system have not been reported. Here we will show how this transition affects the tumor- 
growth and how the effect depends on the initial state of tumor. Our results are inconsistent 
with the zero-dimensional reports that suggest the fluctuation of growth rate always puts 
the tumor at a disadvantage [10, 11]. 

The tumor-growth under immune surveillance can be described by means of insect out- 
break model [1, 16, 17], which in non-dimensional units is given by 



where u is the population of tumor cells; r is their linear per capita birth rate and K is the 
carrying capacity of the environment, respectively. (5v?'/{l -\- v?) quantifies the abilities of 
immune cells to recognize and attack tumor cells. In general, chemotherapy can lead to a 
fluctuation of tumor growth, simply a fluctuation of tumor growth rate r. If considering the 
space of tumor-growth, the growth rate r in Eq.(l) should be rewritten as ro + ^i(^), where 
^j(i) is the Gaussian noises, white in time and space, with zero mean and autocorrelation 
deflned by — 0, — '2a'^6ijS{t — t'), in which cr^ is the noise level and i,j 

are lattice sites. The equivalent stochastic differential equation of Eq.(l) will be. 



here n{i) is the set of the 2d nearest neighbors of site i, d and D are the spatial dimension 
and the diffusion coefficient, respectively. 

Equations of this kind are general and cover different tumor growth and diffusion phe- 
nomena, especially nonequilibrium growth. Wc would like to track down the existence of 
nonequilibrium phase transition induced by multiplicative noise, in systems described by 
these equations. Such a phase transition is characterized by the appearance of multiple 
steady state probability distributions Psti{ui}), which has been applied successfully in nu- 
merous stochastic problems [18, 19]. If set f{ui) — roUi{l — Ui/K) — I3u1/{1 -\- uf), and 
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g{ui) — Ui{l — Ui/K), one will obtain the equivalent Fokker-Planck equation of Eq.(2), 

dp{{u^},t) _ ^[A{uM{u^},t)] ^'[B{uM{u^},t)] . . 

dt dui dul ^ ' 

in which 

A{ui) = f{ui) + a^g{ui)g'{ui) + ^ X] (wi - uj) 

jen{i) 

B{u,) = a^g\u,) (4) 

For simplicity of notation, we drop the subscript i. The stationary solution to Eq.(3) is 
given to be 

r2 r. f{u)-^4g{u)g{u)-D[v-E{v)\ ^ 
Pst{u) = Z exp[— dv J (5) 



Jo 

where Z is a normalization constant, and 

E{v) = {vi\vj) = j VjPst{vj\vi)dvj, (6) 

represents the steady state conditional average of vj at neighboring sites j e n{i), given the 
value Vi at site i. 

Using the Weiss mean field approximation [20, 21], neglecting the fluctuation in the neigh- 
boring sites, i.e., E{v) — (u), independent of v, and imposing the self-consistent requirement 
m — (u), we obtain 



J^"" upst{u,m)du 
Jq°° Pst{u,m)du 
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FIG. 1: m as a function of given by Eq.(7). The points correspond to the intersection of curves 
in Fig.l. The critical immune coefficients are = 2.156 and (3^2 = 2.209, respectively, which 
divide the state of a tumor into three levels: excited (E), sub-excited (S) and non-excited (N). 



The solution, m, of the self-consistency equation is the intersection point between F{m) = m and F{m) 

noise level cr^ = 8.0 x 10^. 

The numerical solution of this last equation for parameter values ro = 1.0, D — 0.01, 
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FIG. 2: Stationary probability distributions of average population of tumor cells for different noise 
intensities and immune coefficients. The parameters are (a) (3 = 2.12, = 0.01, (b) (3 = 2.30, cr^ = 
0.01, (c) P = 2.12, cr^ = 0.40, (d) (3 = 2.30, = 0.40. 

and (j^ = 8.0 x 10~^ is shown in Figs.l and 2. The solution, m, as a function of immune 
coefficient, f3, is obtained by the intersection point between F{m) = m and F{m) = y{m) 
(here y{m) represents the function in the middle position of Eq.(7)). Obviously, the average 
populations of tumor cells exhibit monostable state for low and high values of (3, but un- 
stable state for intermediate value of (3. The critical points are /S^^ = 2.156 and (3^2 — 2.209, 
which divide the states of tumor into three levels: excited state (E), sub-excited state (S), 
and non-excited state (N). Here E and N correspond to stable states but S represents an 
un-stable state, which has two or three possible values. This result means the state of tumor 
is determined by the immune coefficient for low value of noise intensity. 

When the noise level cr^ increases, what will happen? To answer this question, we consider 
E and N, respectively, shown in Fig.3, the stationary probability distributions Pst{u) change 
from monostable state to bistable state with increasing noise intensity and more quantitative 
results are given by Fig.4. For a tumor with excited state, shown in Fig.4, when noise level 
increase, its growth can be hold back to a sub-excited state. Conversely, for the non-excited 
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FIG. 3: m as a function of given by Eq.(7). The points are obtained by a method as same as 
for Fig.2. 

tumor, noise can lead the tumor to sub-excited state or even cancerization. This theoretical 
results are confirmed by corresponding simulations of a one-dimensional system, shown in 
Fig. 5, obtained through a numerical integration of the set of stochastic differential equations 
(2) [22, 23]. In the simulation, we consider three sizes but not find one-dimensional finite 
size effect. It is an important future work to analyze multi-dimensional phase transition of 
tumor system in such a homogeneous circumstance. 

In conclusion, we have found strong evidence for the existence of a noise-induced different 
nonequilibrium phase transitions of tumor growth, in which whether the noise advantage 
the tumor depends on the initial state of tumor. When the tumor is excited, noise induces 
a decay. On the contrary, if the tumor is inactive, the noise can stimulate its growth. 
Provided that the noise results from the treatment as chemotherapy, our results suggest 
that estimating the state of a tumor is a crucial work just before treatment begins. 

This work was partially supported by the National Natural Science Foundation of China 
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FIG. 4: One dimensional simulation for the relationship between m and a^. The parameters are 
same as for Fig.4 

(Grant No. 60471023). 
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